%% dx = (0.125,0.125)
solution = [
0.998845 1.13166 1.28222 1.45284 1.64604 1.86512 2.11304 2.39388 2.7123 
1.13166 1.28213 1.45276 1.64615 1.86496 2.11342 2.39428 2.71244 3.07318 
1.27963 1.44976 1.6428 0 0 0 2.7079 3.06739 3.47521 
1.44134 1.63285 0 0 0 0 0 3.4552 3.91458 
1.61442 1.8284 0 0 0 0 0 3.86917 4.38437 
1.79592 2.03426 0 0 0 0 0 4.30374 4.87551 
1.9796 2.24221 2.53989 0 0 0 4.18613 4.74212 5.37202 
2.15872 2.44497 2.76956 3.13756 3.55497 4.02787 4.56351 5.16985 5.85644 
2.3256 2.63394 2.98365 3.38016 3.82973 4.33911 4.91604 5.56926 6.30881 
];
exact = [
1 1.13315 1.28403 1.45499 1.64872 1.86825 2.117 2.39888 2.71828 
1.13278 1.28361 1.45452 1.64819 1.86764 2.11631 2.3981 2.7174 3.07922 
1.2807 1.45122 1.64445 0 0 0 2.71123 3.07223 3.48129 
1.44235 1.63439 0 0 0 0 0 3.46001 3.92071 
1.61515 1.8302 0 0 0 0 0 3.87453 4.39042 
1.79517 2.03419 0 0 0 0 0 4.30638 4.87977 
1.97712 2.24036 2.53867 0 0 0 4.18555 4.74285 5.37436 
2.15447 2.44133 2.76639 3.13473 3.55212 4.02507 4.56101 5.1683 5.85645 
2.31978 2.62865 2.97865 3.37526 3.82467 4.33391 4.91097 5.56486 6.30581 
];
err = [
-0.00115514 -0.00149001 -0.00180591 -0.00215546 -0.00267867 -0.00312405 -0.00396032 -0.00499141 -0.00598669 
-0.00112365 -0.00147414 -0.00176033 -0.00203607 -0.00268068 -0.00288679 -0.00381551 -0.00495563 -0.00603557 
-0.00107094 -0.00145647 -0.00165033 0 0 0 -0.00333585 -0.0048401 -0.00607955 
-0.00100788 -0.00154696 0 0 0 0 0 -0.00481255 -0.00613385 
-0.000723221 -0.00180223 0 0 0 0 0 -0.00536279 -0.00604999 
0.000749905 6.92555e-05 0 0 0 0 0 -0.00263444 -0.00425489 
0.00248333 0.00184249 0.00122061 0 0 0 0.000575326 -0.000734782 -0.0023347 
0.00424864 0.00364029 0.00317362 0.00282392 0.00285319 0.00279936 0.00250369 0.00154682 -8.30071e-06 
0.00581906 0.00528755 0.0049996 0.00490645 0.00506088 0.00519161 0.00507669 0.00440788 0.00300235 
];
